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ABSTRACT 

We investigate the ability of state-of-the-art redshift-space distortions models for the 
galaxy anisotropic two-point correlation function £,{r^_,r\\ ), to recover precise and unbiased 
estimates of the linear growth rate of structure /, when applied to catalogues of galaxies char- 
acterised by a realistic bias relation. To this aim, we make use of a set of simulated catalogues 
at z = 0.1 and z = 1 with different luminosity thresholds, obtained by populating dark-matter 
haloes from a large N-body simulation using halo occupation prescriptions. We examine the 
most recent developments in redshift-space distortions modelling, which account for non- 
linearities on both small and intermediate scales produced respectively by randomised mo- 
tions in virialised structures and non-linear coupling between the density and velocity fields. 
We consider the possibility of including the linear component of galaxy bias as a free parame- 
ter and directly estimate the growth rate of structure /. Results are compared to those obtained 
using the standard dispersion model, over different ranges of scales. We find that the model of 
Taruya et al. (2010), the most sophisticated one considered in this analysis, provides in gen- 
eral the most unbiased estimates of the growth rate of structure, with systematic errors within 
±4% over a wide range of galaxy populations spanning luminosities between L > L* and 
L > 3L* . Accounting for the scale-dependence of galaxy bias plays a crucial role in recov- 
ering an unbiased estimate of / when fitting quasi non-linear scales. Its impact is found to be 
more severe for highly-biased tracers such as Luminous Red Galaxies, for which systematic 
effects in the modelling might be more difficult to mitigate and have to be further investigated. 
Eventually, we test the impact of neglecting the presence of non-negligible velocity bias with 
respect to mass in the galaxy catalogues. Our tests show that this can produce an additional 
systematic error of the order or lower than 2% on average, comparable to the statistical errors 
we aim at achieving with future high-precision surveys. 
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1 INTRODUCTION 

The structure in the Universe grows through the competing effects 
of universal expansion and gravitational instability. For this reason, 
the large-scale spatial distribution and dynamics of galaxies, which 
follows in some way those of mass, provide fundamental informa- 
tion about the expansion history and the nature of gravity. In gen- 
eral, galaxies are biased tracers of the underlying mass distribution. 
However, they are sensitive to the same gravitational potential and 
their motions keep an imprint of the rate of structure growth. One 
manifestation of this is the observed anisotropy between the clus- 
tering of galaxies along the line-of-sight and that perpendicular to 
it in redshift space. These anisotropies or distortions are caused by 
the line-of-sight component of galaxy peculiar velocities affecting 
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the observed galaxy redshifts from which distances are measured. 
In turn, the large-scale coherent component of galaxy peculiar mo- 
tions is the fingerprint of the growth rate of structure. 

By mapping the large-scale structure over scales which retain 
this primordial information, galaxy spectroscopic surveys have be- 
come one of the most powerful probes of the cosmological model. 
A specific application of spectroscopic surveys involves recover- 
ing cosmological information on the expansion history H{z), by 
measuring the shape of the power spectrum (e.g. Tegmark et al. 
2004; Cole ct al. 2005; Pcrcival et al. 2007; Reid et al. 2010) and 
tracking the baryonic acoustic oscillations (BAO) feature in the 
power spectrum or in the two-point correlation function at differ- 
ent redshifts (e.g Eisenstein et al. 2005; Pcrcival ct al. 2010; Kazin 
et al. 2010; Blake et al. 2011b). However, measurements of H{z) 
alone, either from BAO or Type la supernovae, cannot discrim- 
inate dark energy from modifications of General Relativity (e.g. 
CaiToll et al. 2004), in order to explain the observed recent accel- 
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eration of the expansion of the Universe (Riess et al. 1998; Perl- 
muttcr ct al. 19' ). This degeneracy can be lifted by measuring the 
growth rate at different epochs (Peacock et al. 2006; Albrecht et al. 
2006). Indeed, scenarios with similar expansion history but differ- 
ent gravity or type of dark energy will have a different rate of struc- 
ture growth resulting from different effective gravity strength in 
action. This makes redshift-space distortions measured from large 
spectroscopic surveys a very efficient probe to test cosmology, at 
the same level as BAO and cosmological microwave background 
(CMB) anisotropics. In fact, although this effect is known since the 
late eighties ( iaiser 1987), its usefulness as a probe of dark energy 
and modified gravity has been realised only recently (Guzzo ct al. 
2008). 

Measuring the growth rate of structure from redshift-space 
distortions is however non trivial. The linear theory formalism for 
the power spectrum was first derived by Kaiser (1987) (see Hamil- 
ton 1992, for its configuration-space counterpart). Its validity is 
however limited to very large scales as it lacks a description of 
small-scale non-linear fluctuations. This model has been extended 
to quasi- and non-linear scales in the early nineties using the earlier 
ideas of the "streaming model" (Peebles 1980), in which the lin- 
ear correlation function is convolved along the line-of-sight with a 
pairwise velocity distribution (Fisher et al. 1994; Peacock & Dodds 
1994). This enables one to approximately reproduce the Fingers-of- 
God small-scale elongation ( ickson I97_). Fitting functions cal- 
ibrated on simulations have also been proposed for this purpose 
(Hatton & Cole 1 999; Tinker cl al. 2006). Such extension of the lin- 
ear model, usually referred as the "dispersion model", has been ex- 
tensively used to measure the growth rate of structure / or the dis- 
tortion parameter /3 = from redshift surveys, using measure- 
ments of both redshift-space correlation function (Peacock et al. 
2001; Hawkins et al. 2003; Zehavi et al. 2005; Ross et al. 2007; 
Okumura et al. 2008; Guzzo et al. 2008; da Angela et al. 2008; 
Cabre & Gaztaiiaga 2009a,b; Samushia et al. 2011b) and power 
spectrum (Percival et al. 2004; Tegmark et al. 2004, 2006; Blake 
et al. 2011a). We refer the reader to Hamilton (1998) for a review 
of older studies. Although generally the dispersion model is found 
to be a good fit on linear and quasi-linear scales (Percival & White 
2009; Blake et al. 201 la), it breaks down in the non-linear regime 
(Taruya et al. 2010; Okumura & Jing 2011). In particular, it has 
been shown that it introduces systematic errors of about 10 — 15% 
on the growth rate parameter (e.g. Taruya et al. 2010; Okumura & 
Jing 201 1; Bianehi et al. 2012), of the order of or greater than the 
statistical errors expected from on-going and prospected very large 
spectroscopic surveys such as WiggleZ (Drinkwater et al. 2010), 
GAMA (Driver et al. 201 1), VIPERS (Guzzo et al. 2012), BOSS 
(White et al. 2011), Euclid (Laureijs et al. 2011), or BigBOSS 
(Schlegel et al. 20 1 1 ). In particular, Euclid, the recently selected 
ESA dark energy mission, should be able to constrain the growth 
rate at the percent level (Wang et al. 2010; Samushia et al. 201 I '; 
Majerotto et al. 2012). There is therefore a strong need to go beyond 
the dispersion model, in order to bring systematic errors below this 
expected level of precision, i.e. measure the growth rate parameter 
in an unbiased way. 

Work in this direction started since quite some time, concen- 
trating first on describing the redshift-space clustering and dynam- 
ics of dark matter. Scoccimarro (2004) demonstrated that the dis- 
persion model gives rise to unphysical distributions of pairwise ve- 
locities and proposed an ansatz that accounts, to some extent, for 
the non-linear coupling between the velocity and the density fields. 
This latter model has been shown to provide a better match to the 
observed redshift-space power spectrum in dark matter simulations 



(Scoccimarro 2004; Jennings et al. 20 M ) and has later on been re- 
fined by Taruya et al. (2010). Further approaches have also been 
proposed while completing this paper (Seljak & McDonald 201 1; 
Reid & White 20 1 1 ), but we shall not discuss them in this analysis. 

All mentioned advanced redshift-space distortions models 
have been tested so far only on the redshift-space power spectra 
of dark matter and dark matter haloes, as extracted from large N- 
body simulations (Kwan et al. 2011; Reid & White 2011; Oku- 
umra et al. 201 1; Nishiiriichi & Taruya 201 1). This is quite differ- 
ent from a real survey, in which the most useful tracers of mass, the 
galaxies, are in general biased with respect to the underlying den- 
sity field through a bias which is generally non-linear and scale- 
dependent. The performance of redshift-space distortions models 
applied to galaxy populations with a priori unknown biases, has to 
be precisely investigated. This is the aim of this paper, in which 
we confront non-linear models of redshift-space distortions for the 
anisotropic two-point correlation function in the case of realistic 
galaxy samples. This is done in the framework of the concordant 
ACDM cosmological model. We study the ability of these mod- 
els to recover the linear growth rate of structure and their range 
of applicability. Furthermore we investigate the effects of galaxy 
non-linear spatial and velocity biases and quantify how the latter 
affect the estimated linear growth rate for differently biased galaxy 
populations. 

The paper is organised as follows. In Sect. 2, we present the 
redshift-space distortions formalism and how models can be imple- 
mented in practice. In Sect. 3, we present the comparison between 
the different models and study the effect of galaxy non-linear bias. 
In Sect. 4, we investigate the impact of neglecting galaxy veloc- 
ity bias in the modelling. In Sect. 5, we summarise our results and 
conclude. 



2 REDSHIFT-SPACE DISTORTIONS THEORY 
2.1 Fourier space 

The peculiar velocity v alters objects apparent comoving position 
s from their true comoving position r, as 



aH{a) 



(1) 



where H{a) is the Hubble parameter, a is the scale factor, and ey is 
the line-of-sight unit vector. The redshift-space density field 5" (s) 
can be obtained from the real-space one by requiring mass conser- 
vation, i.e. [1 + S'is)] d^s = [1 + 5{r)] cPr, as 



5\s) = [1^5{r)] 



(Ts 



1. 



(2) 



In the following we shall work in the plane-parallel approximation 
and in this limit, the Jocobian of the real- to redshift-space trans- 
formation can be simply written as. 



(3) 



where we defined U||(r) = —v\\{r)/{faH{a)) with / being the 
linear growth rate. The linear growth rate parameter is defined as 
the logarithmic derivative of the linear growth factor D{a) and 
given by /(a) = dlnD/dlna. To a very good approximation it 
has a generic form (Wang & Steinhardt 1998; Linder 2005), 



/(a) ~ n^iaf 



(4) 
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where 



(5) 



In this parametrisation, while Q.m characterises the mass content in 
the Universe, the exponent 7 directly relates to the theory of gravity 
(e.g. Lindcr 20(! ). General Relativity scenarios have 7 ~ 0.55. 

From Eq. 2 and Eq. 3, one can write the redshift-space density 
field as, 



(6) 



One usually assumes an irrotational velocity field for which 
■it||(r) — d\\lS.~^0{r) and where 9{r) = V ■ v{r) is the veloc- 
ity divergence field and A denotes the Laplacian. In that case Eq. 6 
can be recast, 

5^{s) = {5{r) + fdl^-^e{r)) (l - /a^A-^e(r))"' . (7) 

In Fourier space, it is noticeable that 9y = (fcy/fc)^ — fi^ 
with fj. being the cosine of the angle between the line-of-sight and 
the separation vector. Therefore, one can write the redshift-space 
density field ( >coccimarro et al. 1999) as, 

S%k,^i)= f 7^e-*^5=(s) (8) 



(27r)3 



-ik-r —ikfjj, 



(27r)3 

and the redshift-space power spectrum as. 



[S{x) + fe{x)] (9) 



[&{x) + t?fe[x)mx') + ^^fe{x')]) 



(10) 



where in the latter equation. Any = ^ "11(2;') and r — 

X — x' . The redshift-space power spectrum given in Eq. 10 is al- 
most exact, the only approximation which has been done is to as- 
sume that all object line-of-sight separations are parallel. This ap- 
proximation is valid for samples with pairs covering angles typ- 
ically lower than 10° (Matsubara 2000). Eq. 10 captures all the 
different regimes of distortions. While the terms in square brackets 
describe the squashing effect or "Kaiser effect" which leads to an 
enhancement of clustering on large scales due to the coherent infall 
of mass towards overdensities, the exponential prefactor is respon- 
sible to some extent for the Fingers-of-God effect (FoG, Jackson 
1972) which disperses objects along the line-of-sight due to ran- 
dom motions in virialised structures. .Scoccimarro (2004) proposed 
a simple ansatz for the redshift-space anisotropic power spectrum 
by making the assumption that the exponential prefactor and the 
term involving the density and velocity fields can be separated in 
the ensemble average. In that case Eq. 10 simplifies to, 

(11) 

where Pss, Pse, Pee are respectively the non-linear mass density- 
density, density-velocity divergence, and velocity divergence- 
velocity divergence power spectra and (t„ is the pairwise velocity 
dispersion defined as. 



1 

6tv^ 



■){k)dk. 



(12) 



It is found that this model captures most of the distortion fea- 
tures predicted by N-body simulations (Scoccimarro 2004; Jen- 
nings et al. 201 1 ) although it breaks down in the non-linear regime 
(Percival & White 2009; Taruya et al. 2010). Note that in the linear 



regime where Pgg — Pse = Pee ~ P and in the limit where ka^ 
tends to zero, one recovers the original Kaiser ( 1 987) formula. 



P=(fc,M) = [1 + 2mV + M 



(13) 



derived from linear-order calculations. 

In principle, the exponential prefactor and the term involving 
the density and velocity fields in Eq. 10, which we will refer to as 
the damping and Kaiser terms in the following, cannot be treated 
separately. Additional terms may arise in Eq. 1 1 from the coupling 
between the exponential prefactor and the velocity divergence and 
density fields. Taruya et al. (2010) proposed an improved model 
that takes into account these couplings, adding two correction terms 
Ca and Cs to Scoccimarro (2004)'s formula such as, 

P''ik,ix)^D{kfia,) [P66{k) + 2fi''fPse{k) + fi^fPee{k) 
+CA{k,ii-J) + CB{k,ii;f)], (14) 

whose perturbative expressions are given in their appendix A. In 
the improved model, the exponential prefactor has been replaced 
by an arbitrary functional form D{kfiav) for which is an ef- 
fective pairwise velocity dispersion parameter that can be fitted for. 
Taruya et al. (2010) showed that while adopting a Gaussian or a 
Lorentzian for the damping function and letting (t„ free, one im- 
proves dramatically the fit to the redshift-space power spectrum in 
large dark matter simulations, particularly on translinear scales. 

The function D{kjj,av) damps the power spectra in the Kaiser 
term but also partially mimics the effects of the pairwise velocity 
distribution (PVD) in virialised systems, which translate into the 
FoG seen in the anisotropic power spectrum and correlation func- 
tion on small scales. This is analogous to the phenomenological 
dispersion model proposed in the early nineties (e.g. l-ishcr et al. 
1994; Peacock & Dodds 1994) in which the linear Kaiser model 
in configuration space (Hamilton 1992) is radially convolved with 
a PVD model to reproduce the FoG elongation on small scales, as 
for the early streaming model ( 'ecbles 1980). 

There is however not any simple general functional form for 
the PVD that matches all scales for all types of tracers. The shape 
of the PVD is found to depend on galaxy physical properties and 
halo occupation (Li el al. 2006; Tinker et al. 2006), and its asso- 
ciated pairwise velocity dispersion to vary with scale, in particular 
at small separations (e.g. Hawkins et al. 2003; Cabre & Gaztafiaga 
2009b). It can be shown mathematically that the PVD is in fact not 
a single function but rather an infinite number of PVD correspond- 
ing to different scales and angles between velocities and separation 
vectors ( coccimarro 2004). In practice however, the use of an ex- 
ponential distribution, a Gaussian or other forms with more degrees 
of freedom (e.g. Tang et al. 201 1; Kwan et al. 201 1) shows to be 
very useful to fit the residual small-scale distortions remaining once 
the large-scale Kaiser distortions are accounted for, unless one is 
interested in modelling the very small highly non-linear scales. 



2.2 Conflguration space 

The redshift-space anisotropic two-point correlation function can 
be obtained by Fourier-transforming the anisotropic redshift-space 
power spectrum as, 

^("^'""^ = / (C^^''''^'^^'^) = (15) 



where v — ry/s, 



r?., and Li denote Legendre poly- 



nomials. The correlation function multipole moments (s) are de- 
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fined as, 



(16) 



(17) 



In practice, it can be convenient to write tiie redsliift-space two- 
point correlation function as a convolution between the Fourier 
transform of the damping function D and that of the Kaiser term 
as. 



where ji denotes the spherical Bessel functions and 
9/ + 1 



(27r)3 

= -D(r II, (T„ ) (rx,r||) 



(18) 



In the case of the model of Eq. 11, the three non-null correlation 
function multipole moments related to the Kaiser term are given 

by, 

il'"{s) = 6*W + f\u{r) + f\iee{r), (19) 
e2"{s) = f\ifhr) + fijifhr), (20) 

= /'|er;W, (21) 

where faa ,(,se, ^ee are the Fourier conjugate pairs of Pss, Pse, Pee 
and (Hamilton 1992; Cole et al. 1994), 



^ Jo 



5_3_ 
2 r3 



^x{r')r''^dr' 



7_5_ 
2 rs 



/ S,x{r')r'^dr' . 
Jo 



(22) 



(23) 



The correlation function multipole moments of the Kaiser term in 
the case of the model of Eq. 14 are given in Appendix A. For the 
latter model we restrain ourselves to use only the first three non- 
null multipole moments, as those of orders I = 6 and / = 8 are 
very poorly defined in our simulated galaxy catalogues. 



2.3 From mass to galaxies 

The models derived in the previous section apply in the case of per- 
fectly unbiased tracers of mass. Real galaxies however are biased 
with respect to mass. Galaxy biasing is generally expected to be 
non-linear, scale-dependent, stochastic, and to depend on galaxy 
type, although it is still poorly constrained by observations. On 
large scales in the linear regime, one expects the bias to be a con- 
stant multiplicative factor to the mass density field as Sg = BlS. In 
that case, it is convenient to replace the growth rate / in the models 
with a "effective" growth rate (or distortion parameter) /3 = 
which accounts for the large-scale linear bias fei of the considered 
galaxies. This simple model is valid on large scales where the bias 
asymptotes to a constant value but breaks down on small non-linear 
scales, where bias possibly varies with scale. Recently, Okumura & 
Jmg (20 1 1 ) showed that the scale-dependent behaviour of halo bias 
can strongly affect the recovery of the growth rate. While some ana- 
lytical approaches have been proposed to include bias non-linearity 



in the model (Ocsjacques & Sheth 20 1 0; Matsubara 20 1 ' ), here we 
follow a different route and assume that the scale dependence of 
bias is known. In fact, the latter can be measured to some extent 
from the data themselves in configuration space, once the shape for 
the underlying non-linear mass power spectrum is assumed. Gen- 
eral arguments may suggest that galaxy motions are also biased 
with respect to the mass velocity field, while observations tend to 
indicate that this bias is small (Tinker et al. 2006; Skibba et al. 
20 1 1 ). In this analysis we will neglect the galaxy velocity bias in 
the models but discuss and quantify its impact on the recovery of / 
in Section 4. 



2.4 Constructing the galaxy redshift-space distortion models 

We will use in this analysis different combinations of Kaiser terms, 
damping functions, and bias prescriptions. Although we will work 
in configuration space, we refer to the different models in this sec- 
tion as their Fourier-space counterpart for clarity. All the models 
we consider take the general form. 



Pg{k, fi) = D{kfiay)PKik, n, b) 



(24) 



where. 



exp(— (fe^CT„)^) 
l/(l + (fc^ta„)2) 



and 



PK{k,^l,h) = 

( b^{k)Pss{k) + 2^i''fb{k)Pss{k) + ^i^fPss{k) (mod. A) 
b\k)Pss(k) + 2^i\fb{k)Pse{k) + ii^fPee{k) (mod. B) 



b^{k)Pss{k) + 2fi''fb{k)Pse{k) + fJ.*fPee{k) 
[ +CA{k,ix;f,b) + CB{k,ii-f,b) 



(mod. C) 



fe(fc) = 



bL 



bLbNL{k) 



Hereafter, we will refer as the different Pk models to A, 
B, and C. Model A corresponds to the Kaiser (1987) model with 
the non-linear power spectrum instead of the linear one. It as- 
sumes a linear coupling between the density and velocity fields 
such that S oc 0. Model B is the generalisation proposed by Scoc- 
cimarro (2004) that accounts for the non-linear coupling between 
the density and velocity fields, making explicitly appearing the ve- 
locity divergence auto-power spectrum and density-velocity diver- 
gence cross-power spectrum. Finally, model C is an extension of 
model B that contains the two additional correction terms proposed 
by faruya cl al. (2010) to correctly account for the coupling be- 
tween the Kaiser and damping terms. Besides, we will consider 
two deterministic galaxy biasing prescriptions: a constant linear 
bias b{k) — bL and a general non-linear bias which we define 
as b{k) = (Pgg/Pss)^'''^ (k) = bLbNhik), where Pgg is the 
galaxy power spectrum and bNL{k) is the scale-dependent part 
of the bias that tends to unity at small k. We note that the Gaus- 
sian and Lorentzian damping forms that we will consider here, 
correspond respectively to Gaussian and exponential functions in 
configuration-space. 

The redshift-space distortions models necessitate Pss, Pse, 
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Figure 1. The predicted Pgg, Pgg, Peg by CosmicEmu (Lawrence et al. 
2010) and Jennings et al. (201 1) fitting functions (solid and dotted curves), 
and Closure theory (Taruya et al. 2009) (dashed curves) at 2 = in the 
assumed cosmology. In all cases, P^g and Pgg have been divided by a factor 
of 3 and 10 respectively to improve the clarity of the figure. 



In this analysis we use the Pgs provided by CosmicEmu em- 
ulator (Lawrence et al. 20 1 ( ) and the fitting functions of ennings 
et al. (201 1) to obtain Pgg and Pgg from P55. The latter fitting func- 
tions have an accuracy of 5% to fc = 0.2 for both standard AC DM 
and quintessence dark energy cosmological models. In Fig. 1 and 
2 we compare the Pss, Pgg, Pgg obtained in this way with Closure 
Theory 2-loop analytical predictions at 2; = and z — 1. We find 
that all power spectra agree very well below ~ 0.2 and fc ~ 0.3 
respectively for the two redshifts considered, except in the case of 
Pgg for which they systematically differ by about 10%. While all 
other power spectra match on linear scales, the Pgg fitting formula 
from Jennings cl al. (201 1) stays somewhat below (dotted lines in 
the figures). We find that by multiplying the latter by a factor of 
1.1 one obtains a good match with Closure Theory predictions on 
both linear and non-linear scales (solid lines in the figures). We will 
then adopt this correcting factor in the following when calculating 
the redshift-space distortions models. 

It is noticeable that Closure Theory breaks down at lower fc 
than , ennings cl al. (201 I) fitting functions, indicating that the lat- 
ter are more suitable in practice to describe clustering on the small- 
est scales. It is however important to mention that the validity of 
these fitting functions is limited at large fc (kmax — 0.2, Jen- 
nings et al. 2011). One can therefore reliably describe the mod- 
els involving Pgg and Pgg only down to scales of r ~ Tr/kmax — 
10.5 Mpc. On smaller scales, Pgg and Pgg configuration-space 
counterparts, ^gg and ^gg, will drop rapidly and their contributions 
to the redshift-space distortions models will be underestimated with 
respect to that of ^gg. 




0.01 0.1 1 

fe[hMpc"'] 



Figure 2. Same as Fig. 1 but for z = 1. 



and Pgg real-space power spectra as input. These can be obtained 
analytically using perturbation theory. Although standard perturba- 
tion theory does not describe well the shape of these power spectra 
on intermediate and non-linear scales, improved treatments such as 
Renormalised Perturbation Theory (RPT, Crocce & Scocciniarro 
2006) or Closure Theory (Taruya et al. 2009) have been shown to 
be much more accurate (see Carlson et al. 2009, for a thorough 
comparison). In particular. Closure Theory predictions are found 
to match large N-body simulation real-space power spectra to the 
percent-level up to fc = 0.2 for z > 0.5 (Taruya et al. 2009). 



3 MODEL TESTING 
3.1 Methodology 

To test the redshift-space distortions models presented in sec- 
tion 2.4 and quantify how well the linear growth rate parame- 
ter / can be recovered from the anisotropic two-point correla- 
tion function ^(r-x,r||), we constructed a set of realistic galaxy 
catalogues. We populate the identified friends-of-friends haloes 
in the MultiDark Run 1 (MDRl) dark matter N-body simula- 
tion (Prada et al. 20 II) with galaxies by specifying the Halo 
Occupation Distribution (HOD). MDRl assumes a ACDM cos- 
mology with = f^dm + ^b, Oa, ^b, h, n, ag) = 
(0.27, 0.73, 0.0469, 0.7, 0.95, 0.82) and probe a cubic vol- 
ume of 1 h"^ Gpc^ with a mass resolution of rup = 8.721 x 
10^ Mq. From haloes at snapshots z = 0.1 and z — 1, we 
built galaxy catalogues based on the current most accurate obser- 
vations of the halo occupation available at these redshifts (Zheng 
et al. 2007; Zehavi et al. 201 1; Coupon et al. 201 1) and create three 
luminosity-threshold samples corresponding lo L > L* , L > 2L* , 
and L > 31/*. In these catalogues, the redshift-space displace- 
ments with respect to real space were reproduced using Eq. 1 and 
the galaxy peculiar velocity information. The details of the proce- 
dure used to create the galaxy catalogues are given in Appendix 
B. 

We measured the anisotropic two-point auto-correlation func- 
tion ^{r±, ry) in the redshift-space catalogues using the Landy & 
Szalay (1993) estimator. Because of the large number of galaxies 
and to keep the computational time reasonable, all pair counts have 
been performed using a specifically developed parallel kd-tree code 
following the dual-tree approach (Moore cl al. 200 1 ). 

The best-fitting parameters for the different models have been 



© 2011 RAS, MNRAS 000, 1-15 



6 S. de la Torre & L. Guzzo 



determined by adopting the usual likelihood function, 



-21n£ 



JVp JVp 



(25) 



where A'j, is the number of points in the fit, A is the data-model 
difference vector, and C is the covariance matrix of the data. The 
likelihood is performed on the quantity y — In (l + (,{r±, ry)), 
rather than simply ^{r±,r^\), as to enhance the weight on large 
more linear scales (see Bianchi et al. 2012, for discussion). 

The determination of the covariance matrix is however trou- 
blesome when fitting two-dimensional correlation functions. Hav- 
ing only one realisation of the simulated catalogues at each red- 
shift, we can use to this end internal estimators, such as blockwise 
bootstrap or jackknife resampling (e.g. Norberg et al. 2009). Us- 
ing the latter method, we find that the maximum number of cu- 
bic sub-volumes that can be extracted without underestimating the 
variances on the scales of interest is 64. This poses a problem, since 
in order to have a proper estimation of the eigenvalues of the co- 
variance matrix, this number should be at least equal to the number 
of degrees of freedom, which in our case ranges between 14397 
and 25597 (i.e. 120^ to 160^ data points minus 3 free parameters) 
depending on the scale interval of the fit. As a result, the covariance 
matrix estimated in this way is degenerate and in the end not very 
useful. We note that this is a general problem for redshift-space dis- 
tortions analysis, when one tries to fit the full anisotropic two-point 
correlation function or power spectrum. In principle this could be 
solved by using a very large number of realisations or alternatively, 
theoretically-motivated analytical forms for the covariance matrix. 
In our case, due to the limited size of the simulation, it is impracti- 
cable to define more sub- volumes than degrees of freedom, as this 
would result in underestimating the covariances by using too small 
sub-volumes. We note however that by estimating (,{r± ,r\\) using 
bins of size 0.5h~^ Mpc in both directions in the jackknife resam- 
plings, we are oversampling the functions so the actual number of 
degrees of freedom may be smaller than that associated to the num- 
ber of data points (e.g. I'ishcr cl al. 1994). We are therefore forced 
to perform our fits ignoring the non-diagonal elements of the co- 
variance, i.e. use the variances only. We verified however on a test 
case that the best-fit values of the parameters obtained in this way 
do not differ significantly when using the full covariance matrix 
based on 64 sub-volumes. 

In all cases we define /, a^, and 6_t as free parameters and 
use different scale ranges in the fit by varying r± from r™™ = 
Ih-i Mpc to rT" = 20h-i Mpc and fixing r^"'' = = 
80h~^ Mpc. We limit ourselves to these latter maximum scales 
to avoid being affected by wide-angle effects that may introduce 
additional distortions at larger scales (e.g. Samushia ct al. 201 lb). 
The statistical errors on the model parameters, and in particular on 
/, have been estimated from the la dispersion of their best-fitted 
values among the 64 resamplings. We assume the amplitude of the 
input power spectra as in the models to be known and fix it to that 
of the simulation. The fiducial values of the growth rate are 
given by Eq. 4. 
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Figure 3. Top: relative systematic error on / for L > L* galaxies at z = 1 
in the case of models for which galaxy bias is assumed to be linear. The 
light (dark) shaded band marks the 2% (5%) region around the fiducial 
value. Bottom: corresponding Icr statistical errors on /. 
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Figure 4. Same as Fig. 3 but at redshift z = 0.1. Note the very different 
behaviour of all estimators with respect to the z = I case. 



3.2 Varying the Kaiser and damping terms 

Let us first study the effect of using different combinations of 
Kaiser term and damping function with the assumption that galaxy 
bias is linear and scale-independent. We use the full galaxy cata- 
logues at 2: = 1 and 2 = 0.1 and estimate the statistical and sys- 
tematic errors on the growth rate with the different models, varying 



the minimum perpendicular scale used in the fit, r_L. In this and 
the following section, we consider simulated galaxies with lumi- 
nosities L > L* (see Appendix B for details), having linear bi- 
ases of 6l = 1.34 and b_L ~ 2.01 respectively at z = 0.1 and 
z = 1. These values are realistically close to current observations 
(e.g. Norberg et al. 2001; Polio et al. 2005; Coil et al. 2006; Zehavi 
etal. 2011). 
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3.2. 1 Statistical and systematic errors 

The systematic errors on the growth rate at 2 = 1 are significant us- 
ing any of the model variants, and depends on the chosen minimum 
scale for the fit, r^™, as shown in Fig. 3. Over the explored range 
of r^'". A/// varies between —8% and +6% in the case of mod- 
els A and B. The expected most accurate model, i.e. model C with 
exponential damping (C-EXP), shows a similar trend with a system- 
atic error varying between —6% and +10% and between —21% 
and +9% when using a Gaussian damping (C-GAUSS). We note 
that model A with exponential damping (A-EXP) applied to scales 
r™"' < 5h~^ Mpc, which is one of the most commonly used 
model in the literature, systematically underestimates the value of 
/, in agreement with recent analyses (e.g. Bianchi ct al. 2012). 

When looking at the z = 0.1 case in Fig. 4, model per- 
formances and associated scale-dependence of systematic errors 
change significantly. All models now tend to systematically under- 
estimate the growth rate by about 10% for r™'" > 5h"^ Mpc. The 
notable exception is model C, and more particularly model C-EXP, 
which performs very well and recovers the fiducial value of / with 
less than 3% of systematic error, even including scales as small as 

This is consistent with the power spectrum analysis of Kwan 
et al. (201 I), who show that for dark matter only at z = 0, C- 
GAUSS is the least biased model, with a systematic error of —4% 
when fitting up to kmax ~ 0.1. Our tests show however that 
for galaxies, model C-EXP is less biased than C-GAUSS. In fact 
the choice of damping function has only a significant impact on 
model's ability to handle small scales, with the difference dimin- 
ishing with increasing r^^" given the similar asymptotic behaviour 
of the two functional forms. Conversely, we note that the Gaussian 
damping produces in general slightly lower statistical errors than 
the exponential damping. These tend also to be about 15% smaller 
for models A and B than for model C. 

The important result of this section is the clear difference in 
the systematic error on / produced by the very same models when 
applied to samples at different redshifts. In the next section we shall 
argue that this can be simply understood in terms of the evolution 
(weakening) of the scale-dependence of bias between z = 1 and 
z — 0.1 for L > L* galaxies. 
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Figure 5. Measured § (rx , ) and associated models for L > L* galaxies 
at 2 = 1. In each panel the dotted, dashed, and solid curves correspond re- 
spectively to model A, B, and C with exponential damping and linear bias, 
while the contours correspond to the measured g (rx , rn ) in the galaxy cat- 
alogue. The top panel shows the best-fitting model when (/,(t^,6j^) param- 
eters are allowed to vary, while the bottom panel shows the fiducial predic- 
tion of the models. We note that in the latter case CTi, is fixed to its linear 
value. In this figure, the measured § (rx , f"!! ) is smoothed using a Gaussian 
kernel of size 0.5h~^ Mpc. 



3.2.2 Fidelity in reproducing the anisotropic two-point 
correlation function 

We visually compare in Fig. 5 and Fig. 6 the model correlation 
functions at z = 1 and 2: = 0.1 using (a) the fiducial values of 
the parameters (/, (t„, 6l) and (b) their best-fit values obtained 
with the chosen model, to those measured in our simulated cata- 
logues. For case (a), (t„ is set to linear theory prediction which is 
given by Eq. 12, replacing Peg by the linear power spectrum. We 
limit this comparison to the case with exponential damping and 
rx™ ~ 5h~^ Mpc. The galaxy bias is assumed to be linear and 
scale-independent. 

As shown in the top panels, when all parameters are fixed to 
their fiducial value, model C (solid line) provides a very good de- 
scription of the observed £,(r±,r\\) in the catalogues. Conversely, 
model A and B produce contours that at intermediate separations 
(rx < 40h"^ Mpc, ry < 25 - 30h"^ Mpc) are less squashed 
along the line-of-sight. When (/, (t„, 6_l) are allowed to vary 
(lower panels), all models are generally capable to achieve a good 
fit to £,{r^,r\^\) above rx ~ lOh^^ Mpc. For model A and B, 
this is obtained through a lower damping than predicted by lin- 



ear theory. This allows one to reproduce the significant squashing 
along the line-of-sight seen in the catalogues, but fails to model 
the FoG elongation at smaller r±. It should also be noted that at 
rx < 40h~^ Mpc and ry > 25 - 30h~^ Mpc, model A predicts 
a much higher amplitude of the correlation function than models B 
andC. 

In general, low values of can balance the deficit of power 
on small scales but do not allow to recover the true value of /. 
This is shown in Fig. 7 and 8, where the recovered values of (7„ 
for the different models are compared to linear prediction. One no- 
tices immediately that the trends of the best-fitting values of Ov as 
a function of rx ™ are correlated with those found for / in the dif- 
ferent models: higher values of estimated / imply higher values of 
estimated a^. These results evidence the presence of a significant 
correlation among the parameters, in particular between / and a„. 
In this respect, model C tends to be less degenerate than A and 
B. Moreover, model C is able to recover more realistic values of 
(Ji,, of the order of linear theory predictions. Conversely, model A 
provides best-fitting values strongly deviating from linear expec- 
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Figure 6. Same as Fig. 5 but at 2 = 0.1. 

tations, leading to (t„ as small as 100 — 200 km ■ s"'^. These re- 
sults are consistent with the findings of Faruya ct al. (2010) and 
Nishimichi & Taruya (201 i), who compared model A and C in the 
case of dark matter and dark matter halo catalogues, finding a bet- 
ter agreement of the recovered cr„ with linear theory in the case of 
model C. 

3.3 Effect of galaxy scale-dependent bias 

We now let the galaxy bias vary with scale in the models and study 
the effect on the recovery of the growth rate parameter. In general, 
the galaxy bias in configuration space can be defined as. 



fe(r) 



1/2 



hLbNhir) 



(26) 



where is the galaxy real-space auto-correlation function and 
bNL{r) is the non-hnear scale-dependent part of the bias. It is im- 
portant to stress that £,gg{r) is directly measurable from observa- 
tions by deprojecting the observed projected correlation function 
w{r±_) (Saunders et al. 1992) as, 



1 r°° dw{rj_) 
TV .L dr± 



( 2 



71- 5^ 



-1/2 



dr I 



(27) 
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Figure 7. Recovered values of Ci, for the different models in the case of 
L > L* galaxies at z = 1. The symbol definition is the same as in Fig. 
3. The connected points correspond to the case with scale-dependent bias, 
while others correspond to that with linear bias. The latter points are slightly 
shifted horizontally for clarity. 
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Figure 8. Same as Fig. 7 but at z = 0.1. 



where Wi — wij-i) in the logarithmic bin centred on r±.i. This 
procedure allows one to correctly recover the shape of €,gg{r) up to 
about 30h~^ Mpc (e.g. Saunders ct al. 1992; Cabre & Gaztaiiaga 
2009b). In the following we will therefore make the assumption 
that ^gg(r) is known and use the measured real-space £,gg{r) from 
the simulated catalogues to infer hNhir) to be used in the models. 
In fact, it is not necessary to know the exact shape of (,gg(r) on 
scales larger than about 20 — 30 h^^ Mpc, where one generally 
finds the galaxy bias to be almost scale-independent and can thus 
safely assume bNL{r) = 1. A notable exception is that of more 
non-linear objects, for which the scale dependence may extend to 
larger scales (see section 3.3.1). 

Fig. 9 shows the non-linear scale-dependent component of 
galaxy bias, hNhij'), for the different galaxy populations in our 
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Figure 9. The scale dependence of galaxy bias at z = 0.1 and z = 1.0, 
for the different galaxy populations considered in this work (see inset). It is 
defined as 6jvL(r) = [^gg{r)/ (bl^ssir))]^^'' ■ 
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Figure 11. Same as Fig. 10 but at z = 0.1. 
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Figure 10. Top: the relative systematic error on / for L > L* galaxies at 
2 = 1, when the bias scale dependence as measured from the simulation is 
included. The light (dark) shaded band marks the 2% (5%) region around 
the fiducial value. Bottom: conesponding Icr statistical errors. 



Ih^^ Mpc < r < 20h^^ Mpc, the strength of the effect increas- 
ing for more luminous galaxies. 

Let us come back to our original L > L* catalogues and re- 
peat the analysis of the previous section now including the scale 
dependence of galaxy bias shown in Fig. 9. The new statistical and 
systematic errors on / estimated from our simulated catalogues 
are shown in Figs. 10 and 11. In general, one sees that includ- 
ing the bias scale-dependence information improves the estima- 
tion of /, in particular when using the smallest scales in the fit 
(r^'" < 5h~^ Mpc), as expected. This is particularly true atz — 1 
although at z — 0.1 the effect is almost negligible. This is clearly 
related to the scale dependence of the bias being stronger at earlier 
epochs, as visible in Fig. 9. This is a very important result, as it im- 
plies that neglecting bias scale-dependence in analysis of deep red- 
shift surveys, would produce estimates of / which are differently 
affected by systematic errors at different redshifts. Accounting for 
this dependence on scale tends also to reduce the dependence of the 
systematic error on the minimum fitted scale: the retrieved value 
is almost constant down to r^*" — Ih^^ Mpc for all consid- 
ered models. Finally, it is important to note that model C-EXP pro- 
vides the best performance at both considered redshifts, recovering 
the fiducial value of the growth rate with less than 3% systemat- 
ics while using clustering information down to r± — 5h~^ Mpc. 
Remarkably, including the bias scale-dependence also reduces the 
statistical error on / by about 15% for all models. 



simulated catalogues at the two reference redshifts considered, 
z = 1 and z — 0.1. In the previous section we considered only 
catalogues of galaxies with L > L* , while in this figure we in- 
troduce more extreme galaxy populations, which we analyse in the 
following section. To define bNL(r), the linear bias 6i, has been 
determined for each galaxy population by minimising the differ- 
ence between and h\^ii on scales above r = 10h~^ Mpc. 
It is evident from this figure that non-linearities in the galaxy bias 
produce variations up to 40% in the real-space clustering on scales 



3.3. 1 The case of highly-biased galaxies 

One might naively expect highly biased objects to be more suitable 
to constrain / from observations, by virtue of their higher clus- 
tering signal that should in principle reduce statistical uncertainty 
(see Bianchi ct al. 2012). However, highly biased objects, which 
are more likely to reside in the most massive haloes, have under- 
gone a stronger non-linear evolution, this explaining their stronger 
bias scale-dependence seen in Fig. 9. As such, the inclusion of bML 
in the models may become even more critical for these objects. 
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Figure 12. Top: relative systematic error on / for different galaxy popula- 
tions at different redshifts (see inset). In all cases we use model C-EXP and 
include the scale dependence of bias. The light (dark) shaded band marks 
the 2% (5%) region around the fiducial value. Bottom; corresponding la 
statistical errors. 



if one's goal is to accurately measure the growth rate parameter 
(e.g. Cabre & Gaztanaga 2009a). In this section, we explicitly test 
this hypothesis, extending the model comparison to the regime of 
highly biased galaxy populations. 

Let us first define higher-luminosity galaxies from our sim- 
ulated catalogues. We shall consider two sub-samples, defined as 
having galaxies with L > 2L* and L > 3L*. We complement 
these with a catalogue of Luminous Red Galaxies (LRG) drawn 
from the "LasDamas" suite of simulations, meant to accurately re- 
produce the galaxy clustering in the SDSS-DR7 release (McBride 
et al., in preparation). More precisely, we shall use 100 mock real- 
isations of what are defined as "faint LRG" (Air < —21.2). 

We compare in Fig. 12 the relative systematic error on the 
growth rate obtained for these highly biased galaxies. We use here 
only model C-EXP and include the scale-dependence of bias as de- 
scribed in Section 3.3. We find that the trends are similar to those 
observed for L > L* galaxies, with a dependence on r™'" that 
varies with the considered redshift. Overall, the systematic error 
remains within ±5% for rT"" = 5 - lOh"^ Mpc for all con- 
sidered galaxy populations, but the LRG sample. In this case the 
growth rate is overestimated by on average about 10%, suggest- 
ing additional non-linear effects that are not accounted for in our 
models. A cautionary remark in interpreting this result is that the 
simulated LRG sample is the only one which has not been built by 
us on purpose for this work. Although this sample has been already 
used for other investigations, we have no way to verify the details 
of the HOD implementation and its impact on our results. Assum- 
ing the effect is indeed real, one may hope however to reduce the 
systematic error by including in the fit scales greater than our cur- 
rent upper limit of r±_ = 80h~^ Mpc (Samushia et al. 201 la), at 
the expense of a worse precision in the measurement. 

These results suggest that for highly biased galaxies, addi- 
tional systematic effects arise which could be due to an incorrect 



inclusion of scale-dependence of bias or a too simplistic approach 
for the FoG. From a theoretical point-of-view, bias non-linearity 
changes the correction terms Ca and Cb in model C, which might 
need to be modified to properly include galaxy bias non-linearities 
and scale-dependence (see Tang et al. 201 1; Nishimichi & Taruya 
201 1). We plan to investigate in more details these aspects in a 
future paper. Finally, we note that FoG modelling could be poten- 
tially improved by adding more freedom in the damping function 
(e.g. Rwan ct al. 201 1) or including scale dependence in the pair- 
wise velocity dispersion (e.g. Hawkins et al. 200.i), but at the price 
of increasing the statistical error on the growth rate estimate. 



4 EFFECT OF GALAXY VELOCITY BIAS 

In the framework of understanding the impact of non-linear effects 
on the accuracy of growth rate estimates from redshift-space distor- 
tions, we investigate in this section the possible impact of velocity 
bias, an effect which is is usually neglected. The galaxy catalogues 
that we used so far, assume that the radial distribution of satellite 
galaxies inside dark-matter halos follows that of mass, as described 
by a Navarro ct al. (1996) (NFW) radial density profile. Moreover, 
central galaxies have been defined as being at rest at the centre of 
their dark matter halo, inheriting its mean velocity. These assump- 
tions make galaxy velocities unbiased with respect to the mass ve- 
locity field. However, there are some observational evidences that 
galaxies does not exactly follow the same radial distribution as dark 
matter and exhibit some velocity bias. In particular, recent small- 
scale clustering measurements in the SDSS suggest that relatively 
luminous galaxies have a steeper radial density profile than NFW, 
with inner slopes close to —2 and lower concentration parameters 
(Watson et al. 201 1 ). In addition, galaxy groups and galaxy clusters 
analysis tend to indicate that central galaxies might in general not 
be exactly at rest at the centre of the potential well of dark matter 
haloes (e.g. van den Bosch et al. 2005; Skibba ct al. 201 1). Differ- 
ences between the radial distribution of galaxies and that of mass 
implies the presence of additional spatial and velocity biases. This 
has a direct impact on the description of the observed redshift-space 
distortions. In this section we provide a first quantitative assess- 
ment of the systematic uncertainty on / that it can introduce if not 
accounted for in the models. 

To this end, we include in the catalogues some amount of ve- 
locity bias coming from either central or satellite galaxies. For cen- 
tral galaxies we follow van den Bosch et al. (2005) and assign halo- 
centric positions assuming a radial density distribution of the form, 



i(r-|m) oc 



frTv (m) 



r(r + frrv(m)) 



(28) 



where m is the halo mass, is the halo virial radius, and fr is a 
free parameter that controls the amount of spatial and velocity bias 
introduced. This effectively offsets central galaxies from their halo 
centre of mass. By solving Jeans equation one obtains the associ- 
ated one-dimensional velocity dispersion. 
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where 4^{r) is gravitational potential. One can thus define the ve- 
locity bias of central galaxies as. 
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where the halo-averaged velocity dispersions are given by, 
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and the NFW radial density profile is defined as, 
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where Cdm is the dark matter concentration parameter (see Ap- 
pendix B for its definition). Galaxy halo-centric velocities are 
drawn from a Gaussian distribution in each dimension with veloc- 
ity dispersion given by Eq. 29. We choose the value of fr in order 
to have bj*^" (m) = 0.2. The latter is an average value motivated by 
recent observations (e.g. Coziol et al. 2009; Skibba et al. 201 1). 

In the case of satellite galaxies we followed a similar method- 
ology. Based on the recent results of Watson et al, (201 1) in the 
SDSS for Mr < —20.5 galaxies, we radially distribute satellite 
galaxies as to reproduce a generalised radial density profile of the 
form, 
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with Csat {m) = "''"'^^'"^ and 7 = 2. As for central galaxies, we as- 
sign satellite galaxy velocities from the associated one-dimensional 
velocity dispersion. 



at{r\m) = 



1 



P3at{r\m) 
This leads to a velocity bias. 
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where 
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An 



Psatir\m)asatir\m)r dr. (37) 



In that case, (m) slowly increases from ftj"' — 1 at m = 
10^ Mq to fer* 1-2 at m = 10^*^ h'^ Mq. Assuming that 
the haloes are spherical and isotropic, all previous integrals can be 
solved analytically (see e.g. Appendix B and v an den Bosch et al. 
2005). 

Fig. 13 and 14 show the percent variation on the estimated /, 
when a velocity bias is included in the simulated catalogues ac- 
cording to the previously described procedure. The curves show 
the systematic effect of velocity bias coming from either central 
or satellite galaxies, when / is estimated from models B-EXP and 
C-EXP assuming a scale-dependent spatial bias. We find that the 
impact of central galaxy velocity bias is smaller than that of satel- 
lite galaxies: at 2 = 0.1 a velocity bias of bj"^" = 0.2 introduces a 
negative systematic error of about 1%, independently of the model. 
At a = 1 the effect is larger, with a systematic error that can reach 
about -3% over the range Ih"^ Mpc < r™'" < 20h~^ Mpc. 
We should note that the test at z = 1 may be pessimistic, as the 
amount of velocity bias introduced in the catalogues corresponds 
to observational constraints from the local Universe and it is plau- 
sible that at z = 1 the actual velocity bias is lower. In the case of 
satellite galaxies, the introduction of a velocity bias has a stronger 
impact: at z = 1 the systematic error remains within 1 — 2%, while 
at z = 0.1 it is of about 2 — 3%, depending on the model used. In 
the latter case, we note that model C tend to be less affected. 
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Figure 13. Percent systematic variation of the measured growth rate (/"'' — 
jnovb'^ I jnovb induced by including some amount of velocity bias from 
central galaxies in the simulated catalogues. Galaxies with L > L* at the 
two reference redshifts are considered, applying models B-EXP and C-EXP 
with scale-dependent bias. 
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Figure 14. Same as Fig. 13 but when including velocity bias coming from 
satellite galaxies in the simulated catalogues. 



These admittedly simple tests suggest that velocity bias, if not 
accounted for in the models, can introduce additional systematic 
errors of the order of 1 — 3%, i.e. of the same order of statistical 
errors expected to be reachable by future large surveys. Although 
not dramatic, this additional source of systematic error needs to be 
kept in mind and possibly accounted for in future models. In prin- 
ciple, galaxy velocity bias can be included in the models by adding 
an effective velocity bias factor in front of the terms involving the 
velocity divergence field. At first approximation, one could assume 
this velocity bias factor to be constant and set it as a free parameter 
while fitting redshift-space distortions. Introducing more degrees 
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of freedom in the models would however inevitably increase statis- 
tical errors. 



5 SUMMARY AND CONCLUSIONS 

Measurements of the growth rate of structure from redshift-space 
distortions have come to be considered as one of the most promis- 
ing probes for future massive redshift surveys that aim at solving 
the dilemma of cosmic acceleration. A notable example is the sur- 
vey planned for the approved ESA Euclid space mission (Laureijs 
et al. 201 1 ). In this paper, we have investigated in some detail how 
well non-linear effects can be accounted for by current models 
when performing such measurements on galaxy catalogues. The 
question is how to optimally use observations and models, in or- 
der to minimise statistical and systematic errors alike. The actual 
extraction of the redshift-space distortions signal from real data is 
subject to two competing requirements. On one hand, one would 
like to use the simple linear description of redshift-space distor- 
tions induced by large-scale coherent motions, thus limiting the 
measurements and modelling to very large scales. On such scales 
fluctuations are close to be linear and systematic effects appear to 
be reduced (e.g. Samushia ct al. 201 la). On the other hand, the 
clustering signal on those scales is so weak that statistical errors 
on the measured growth rate remain large, even when using sam- 
ples probing large volumes. One may therefore prefer to extend the 
modelling to smaller scales, enlarging the range of analysed scales 
and thus reducing the statistical error. 

The most effective compromise to exploit the size and statis- 
tics of future surveys seems thereby that of including intermediate 
quasi-nonlinear scales, at the expense of a more complicated mod- 
elling effort. As discussed in the introduction, most recent devel- 
opments in this direction have so far concentrated on improving 
the description of non-linear effects in the redshift-space clustering 
of matter fluctuations in Fourier space. Here we have investigated 
how these models perform when applied to catalogues of realis- 
tic galaxies in configuration space, in both the local Universe at 
2 = 0.1 and the more distant Universe at z = 1. For this pur- 
pose, we have reformulated in terms of the anisotropic two-point 
correlation function of galaxies ^(rx,r||), the analytical model 
for the redshift-space anisotropic power spectrum by Scoccimarro 
(2004) as well as the recent improvements proposed by Taruya et al. 
(2010), in addition to the standard dispersion model. In these mod- 
els, we have included the possibility of using a realistic non-linear 
scale-dependent galaxy bias, the latter being in principle measur- 
able from the observations. At variance with the usual habit of fit- 
ting for the distortion parameter /3 = f /bL, v^e have considered 
the possibility of including the linear component of galaxy bias as 
a free parameter and directly estimate the growth rate of structure 
/. Our key findings and results can be summarised as follows: 

(i) When applied to the galaxy anisotropic correlation function, 
Taruya et al. (2010)'s model, the most sophisticated model consid- 
ered in this analysis, generally provide the most unbiased estimates 
of the growth rate of structure /. 

(ii) The inclusion of the scale-dependence of bias in the models 
is an important ingredient for minimising the systematic error on 
/. This is particularly significant in the case of galaxy populations 
with non-negligible non-linearity of the bias, as seen in the z — 1 
samples of the present analysis. 

(iii) Systematic errors vary with the degree of non-linearity in 
the bias of the considered galaxy population, which in turn is a 



function of redshift. Accounting for it is therefore particularly rele- 
vant when "slicing" deep surveys to measure f{z) over a significant 
redshift range. 

(iv) Taruya ct al. (2010)'s model with scale-dependent bias al- 
lows to reach the 4% accuracy on the recovery of the growth rate 
of structure independently of how biased galaxies are. The only ex- 
ception is that of very luminous and biased galaxies such as LRGs, 
for which / can be over-estimated by more than 10%. This shows 
the need to further improving the way bias is included in redshift- 
space distortions models. 

(v) Galaxy velocity bias could represent, if not accounted for, 
an additional source of systematic error. By implementing realistic 
prescriptions for galaxy velocity bias coming from either central 
or satellite galaxies in our simulated samples, we estimate that ne- 
glecting it yields to an underestimate of the recovered / by 1 — 3%. 

Overall, these results emphasize the need for careful mod- 
elling non-linear effects, if redshift-space distortions have to be 
used as a precision cosmology probe. They nevertheless indicate 
promising venues along which to develop further methods to over- 
come systematic biases and support the hope that redshift-space 
distortions in future surveys will provide us with a unique test of 
the cosmological model. 



ACKNOWLEDGMENTS 

We thank Atsushi Taruya for providing us with his codes to com- 
pute Closure Theory predictions and for giving us useful com- 
ments on the manuscript. We also thank Davide Bianchi, John 
Peacock, and Shaun Cole for discussions and suggestions. Finan- 
cial support through PRIN-INAF 2007 and 2008 grants and ASI 
COFISAVP31 10 1/026/07/0 is gratefully acknowledged. 



REFERENCES 

Albrecht, A., Bernstein, G., Cahn, R., et al. 2006, ArXiv Astro- 
physics e-prints 

Bhattacharya, S., Heitmann, K., White, M., et al. 201 1, ApJ, 732, 
122 

Bianchi, D., Guzzo, L., Branchini, E., et al. 2012, in preparation 
Blake, C, Brough, S., CoUess, M., et al. 2011a, MNRAS, 415, 
2876 

Blake, C, Davis, T, Poole, G. B., et al. 2011b, MNRAS, 415, 
2892 

Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, ApJ, 
592,819 

Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, MNRAS, 321, 
559 

Cabre, A. & Gaztaiiaga, E. 2009a, MNRAS, 393, 1183 
Cabre, A. & Gaztaiiaga, E. 2009b, MNRAS, 396, 1119 
Carlson, J., White, M., & Padmanabhan, N. 2009, Phys. Rev. D, 
80,043531 

Carroll, S. M., Duvvuri, V., Trodden, M., & Turner, M. S. 2004, 

Phys. Rev. D, 70, 043528 
Coil, A. L., Newman, J. A., Cooper, M. C, et al. 2006, ApJ, 644, 

671 

Cole, S., Fisher, K. B., & Weinberg, D. H. 1994, MNRAS, 267, 
785 

Cole, S., Percival, W. J., Peacock, J. A., et al. 2005, MNRAS, 362, 
505 

Cooray, A. & Sheth, R. 2002, Phys. Rep., 372, 1 



© 2011 RAS, MNRAS 000, 1-15 



Modelling non-linear redshift-space distortions 13 



Coupon, J., Kilbinger, M., McCracken, H. J., et al. 2011, ArXiv 
e-prints 

Coziol, R., Andemach, H., Caretta, C. A., Alamo-Martmez, K. A., 

& Tago, E. 2009, AJ, 137, 4795 
Crocce, M. & Scoccimarro, R. 2006, Phys. Rev. D, 73, 063519 
da Angela, J., Shanks, T, Croom, S. M., et al. 2008, MNRAS, 

383, 565 

Desjacques, V. & Sheth, R. K. 2010, Phys. Rev. D, 81, 023526 
Drinkwater, M. J., Jurek, R. J., Blake, C, et al. 2010, MNRAS, 
401, 1429 

Driver, S. P, Hill, D. T., Kelvin, L. S., et al. 2011, MNRAS, 413, 
971 

Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 633, 
560 

Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., & Huchra, J. P. 

1994, MNRAS, 267, 927 
Guzzo, L., Pierleoni, M., Meneux, B., et al. 2008, Nature, 451, 

541 

Guzzo, L. et al. 2012, in preparation 
Hamilton, A. J. S. 1992, ApJ, 385, L5 

Hamilton, A. J. S. 1998, in Astrophysics and Space Science Li- 
brary, Vol. 231, The Evolving Universe, ed. D. Hamilton, 185 
Hatton, S. & Cole, S. 1999, MNRAS, 310, 1137 
Hawkins, E., Maddox, S., Cole, S., et al. 2003, MNRAS, 346, 78 
Ilbert, O., Tresse, L., Zucca, E., et al. 2005, A&A, 439, 863 
Jackson, J. C. 1972, MNRAS, 156, IP 

Jennings, E., Baugh, C. M., & Pascoli, S. 2011, MNRAS, 410, 
2081 

Kaiser, N. 1987, MNRAS, 227, 1 

Kazin, E. A., Blanton, M. R., Scoccimarro, R., et al. 2010, ApJ, 
710, 1444 

Kwan, J., Lewis, G. P., & Linder, E. V. 201 1, ArXiv e-prints 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 
Laureijs, R., Amiaux, J., Arduini, S., et al. 201 1, ArXiv e-prints 
Lawrence, E., Heitmann, K., White, M., et al. 2010, ApJ, 713, 
1322 

Li, C, Kauffmann, G., Jing, Y. P, et al. 2006, MNRAS, 368, 21 

Linder, E. V. 2004, Phys. Rev. D, 70, 02351 1 

Linder, E. V. 2005, Phys. Rev. D, 72, 043529 

Majerotto, E., Guzzo, L., Samushia, L., et al. 2012, in preparation 

Matsubara, T. 2000, ApJ, 535, 1 

Matsubara, T. 201 1, Phys. Rev. D, 83, 083518 

Moore, A. W., Connolly, A. J., Genovese, C, et al. 2001, in Min- 
ing the Sky, ed. A. J. Banday, S. Zaroubi, & M. Bartelmarm, 
71-1- 

Navarro, J. R, Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 
563 

Nishimichi, T. & Taruya, A. 201 1, Phys. Rev. D, 84, 043526 
Norberg, P, Baugh, C. M., Gaztaiiaga, E., & Croton, D. J. 2009, 

MNRAS, 396, 19 
Norberg, P, Baugh, C. M., Hawkins, E., et al. 2001, MNRAS, 

328, 64 

Okumura, T. & Jing, Y. P 201 1, ApJ, 726, 5 
Okumura, T., Matsubara, T, Eisenstein, D. J., et al. 2008, ApJ, 
676, 889 

Okumura, T, Seljak, U., McDonald, P, & Desjacques, V. 2011, 
ArXiv e-prints 

Peacock, J. A., Cole, S., Norberg, P, et al. 2001, Nature, 410, 169 
Peacock, J. A. & Dodds, S. J. 1994, MNRAS, 267, 1020 
Peacock, J. A., Schneider, P, Efstathiou, G., et al. 2006, ESA- 
ESO Working Group on "Fundamental Cosmology", Tech. rep. 



Peebles, P. J. E. 1980, The large-scale structure of the universe 

(Princeton University Press, 1980. 435 p.) 
Percival, W. J., Burkey, D., Heavens, A., et al. 2004, MNRAS, 

353, 1201 

Percival, W. J., Nichol, R. C, Eisenstein, D. J., et al. 2007, ApJ, 
657, 645 

Percival, W. J., Reid, B. A., Eisenstein, D. J., et al. 2010, MNRAS, 
401,2148 

Percival, W. J. & White, M. 2009, MNRAS, 393, 297 
Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 
565 

Polio, A., Meneux, B., Guzzo, L., et al. 2005, A&A, 439, 887 
Prada, F, Klypin, A., Yepes, G., Nuza, S. E., & Gottloeber, S. 

201 1, ArXiv e-prints 
Reid, B. A., Percival, W. J., Eisenstein, D. J., et al. 2010, MNRAS, 

404, 60 

Reid, B. A. &. White, M. 2011, MNRAS, 417, 1913 
Riess, A. G., Filippenko, A. V., Challis, P, et al. 1998, AJ, 116, 
1009 

Ross, N. P, da Angela, J., Shanks, T, et al. 2007, MNRAS, 381, 
573 

Samushia, L., Percival, W. J., Guzzo, L., et al. 2011a, MNRAS, 
410, 1993 

Samushia, L., Percival, W. J., & Raccanelli, A. 2011b, ArXiv e- 
prints 

Saunders, W, Rowan-Robinson, M., & Lawrence, A. 1992, MN- 
RAS, 258, 134 

Schlegel, D., Abdalla, F, Abraham, T, et al. 201 1, ArXiv e-prints 

Scoccimarro, R. 2004, Phys. Rev. D, 70, 083007 

Scoccimarro, R., Couchman, H. M. P., & Frieman, J. A. 1999, 
ApJ, 517, 531 

Seljak, U. & McDonald, P 201 1, JCAP, 1 1, 39 

Skibba, R. A., van den Bosch, F C, Yang, X., et al. 2011, MN- 
RAS, 410, 417 

Tang, J., Kayo, I., & Takada, M. 201 1, MNRAS, 416, 2291 

Taruya, A., Nishimichi, T, & Saito, S. 2010, Phys. Rev. D, 82, 
063522 

Taruya, A., Nishimichi, T, Saito, S., & Hiramatsu, T. 2009, 

Phys. Rev. D, 80, 123503 
Tegmark, M., Blanton, M. R., Strauss, M. A., et al. 2004, ApJ, 

606, 702 

Tegmark, M., Eisenstein, D. J., Strauss, M. A., et al. 2006, 

Phys. Rev. D, 74, 123507 
Tinker, J. L., Weinberg, D. H., & Zheng, Z. 2006, MNRAS, 368, 

85 

Toyoda, M. & Ozaki, T. 2010, Computer Physics Communica- 
tions, 181, 277 

van den Bosch, F. C, Norberg, P, Mo, H. J., & Yang, X. 2004, 

MNRAS, 352, 1302 
van den Bosch, F. C, Weinmann, S. M., Yang, X., et al. 2005, 

MNRAS, 361, 1203 
Wang, L. & Steinhardt, P J. 1998, ApJ, 508, 483 
Wang, Y, Percival, W., Cimatti, A., et al. 2010, MNRAS, 409, 

737 

Watson, D. F, Berlind, A. A., McBride, C. K., Hogg, D. W., & 

Jiang, T. 201 1, ArXiv e-prints 
White, M., Blanton, M., Bolton, A., et al. 201 1, ApJ, 728, 126 
Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 201 1, ApJ, 736, 59 
Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2005, ApJ, 630, 1 
Zheng, Z., Berlind, A. A., Weinberg, D. H., et al. 2005, ApJ, 633, 

791 

Zheng, Z., Coil, A. L., & Zehavi, I. 2007, ApJ, 667, 760 



©2011 RAS, MNRAS 000, 1-15 



14 S. de la Torre & L. Guzzo 



APPENDIX A: REDSHIFT-SPACE ANISOTROPIC 
TWO-POINT CORRELATION FUNCTION FOR THE 
TARUYA, NISHIMICHI & SAITO (2010) MODEL 

The redshift-space anisotropic two-point correlation function is 
obtainable by Fourier-transforming the anisotropic redshift-space 
power spectrum P" (k, fi) as, 



where ly' = t\\/ s, r±_ = ^ — , and L; denote Legendre poly- 
nomials. The correlation function multipole moments (s) are de- 
fined as, 



27r2 



k^Pnk)ji(ks), 



(A2) 



where ji denotes the spherical Bessel functions and 



Pnk) = -^ d^,p\k,^l)Ll{^,). (A3) 



In the case of biased tracers of mass, Taruya et al. (2010) 
model for the redshift-space anisotropic power spectrum can be 
written as, 

P=(fe,M) = D{ktia^) [b^Pssik) + 2bfi^fPse{k) 

+tJ,\fPgeik) + CA{k, /, 6) + Ceik, fi; f, b)] 



(A4) 



where b is the spatial bias of the considered tracers and, 
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where functions Amn{r, x), Amnir, x), amn{r, x), Bab{r, x) are 
given in Appendix A of faruya cl al. (2010), P{k) is the linear 
mass power spectrum, Pi2{k) = Psg{k), and P22{k) — Peeik). 
By using the Kaiser term in Eq. A4 (i.e. Eq. A4 without the damp- 
ing function D{kfj,av)) into Eq. A3 and A2, one obtains the cor- 
responding correlation function multipole moments. The non-null 



multipole moments are then given by, 
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where £,Amn and £,Bnab are the Fourier conjugate pairs of PAmn 
and PBnab in Eqs. A5 and A6, and are the correlation function 
multipole moments associated with Px as defined in Eq. A2. For 
orders i = 2, / = 4, Z = 6 and I — 8, the latter can be conveniently 
rewritten as. 



Px as defined in Eq. A2. For 
the latter can be conveniently 

d''M=£x(r)-3X2(r) 
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These identities, which orders I — 2 and I = 4 were already found 
by Hamilton ( 1 992) and Cole et al. ( 1 994), are obtained by using re- 



where. 



X„(r) =r-("+^) f ix{r')r'"dr'. 
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currence relations and integral forms of spherical Bessel functions 
(see Toyoda & Ozaki 20 IC, for details). 

Although the model predicts non-null multipole moments of 
orders I — 6 and I = 8, we neglected these terms in this analy- 
sis, including only correlation function multipole moments ^□(■s), 
^2{s), and (see main text). 



APPENDIX B: HOD GALAXY CATALOGUE 
CONSTRUCTION 

We describe in this appendix the method that we used to create real- 
istic galaxy catalogues from a large N-body dark matter simulation, 
based on the Halo Occupation Distribution (HOD) formalism (e.g. 
Cooray & Sheth 2002). 

We used the dark matter haloes identified using a friends-of- 
friends algorithm with linking length of I = 0. 17 in the MDRl dark 
matter simulation by Prada et al. (2011). Halo masses were esti- 
mated from the sum of dark matter particle masses inside each halo 
after correcting for finite force and mass resolution (Bhattacharya 
et al. 201 ] , Eq. 4). We populated haloes according to their mass by 
specifying the galaxy halo occupation which we parametrise as. 



{Ngai\m) = {Ncen,\m) (1 + {Nsat\m)) 



(Bl) 



where {Ncm\Tn) and {NBat\m) are the average number of cen- 
tral and satellite galaxies in a halo of mass m. Central and satellite 
galaxy occupations are defined as (Zheng et al. 2005), 



{Ncen\m) = 
{Nsat\m) 



1 + erf 

m — Mo 
Ml 



log m - log Mn 



(B2) 
(B3) 



where Mmin, o"iog m, Mo, Ml, and a are HOD parameters. 

We positioned central galaxies at halo centres with probability 
given by a Bernoulli distribution function with mean taken from Eq. 
B2 and assigned host halo mean velocities to them. The number of 
satellite galaxies per halo is set to follow a Poisson distribution with 
mean given by Eq. B3. We assumed that satellite galaxies follow the 
spatial and velocity distribution of mass and randomly distributed 
their halo-centric radial position as to reproduce a Navarro et al. 
(199(i) (NEW) radial profile. 



pNFw{r\m) oc 



Cdv 
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where Cdm is the concentration parameter and (m) is the virial 
radius defined as. 



r„(m) 



3m 



1/3 



(B5) 



In this equation, p{z) is the mean matter density at redshift z and 
Anl = 200 is the critical overdensity for virialisation in our def- 
inition. Consistently with the works of z^chavi cl al, (201 1), Zheng 
et al. (2007), and Coupon et al. (201 1) that we used to infer the 
mean galaxy occupation, we assumed the mass-concentration rela- 
tion of Bullock ct al. (2001), 



Cdm{m, z) = 



Co 



m 



(B6) 



where co = 11, /3 = —0.13, and m, is the non-linear mass scale 
at 2 = defined such as cr(m*,0) = Sc. Here 5c and a{m,0) 
are respectively the critical overdensity (we fixed Sc — 1.686) and 



the standard deviation of mass fluctuations at z = 0. The latter is 
defined as, 

dk k''^P{k,z) 



a (m, z) — 



\W{kR)\ 



(B7) 



where R = [3m/ (47rp(2:))]^''^, P{k, z) is the linear mass power 
spectrum at redshift z in the adopted cosmology, and W{x) is the 
Eourier transform of a top-hat filter 

In order to assign satellite galaxy velocities, we assumed halo 
isotropy and sphericity, and drawn velocities from Gaussian dis- 
tribution functions along each Cartesian dimension with velocity 
dispersion given by (\ an den Bosch et al. 2004), 



1 



pNFw{r\m 
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dr 



1 + ) i(^r/rs) 
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I{x) = 



(BIO) 



where i/'(r) is the gravitational potential, G is the gravitational con- 
stant, fix) — ln(l + x) — x/{l + x), and 

f(t)dt 

We generated catalogues corresponding to different popula- 
tions of galaxies defined such that their luminosity is greater than 
multiples of the characteristic luminosity at both z = 0.1 and 
z — 1. For this purpose, we used the HOD parameters obtained 
from the SDSS survey at 2 ~ 0.1 by /ehavi ct al. (201 i) and 
interpolated the parameter dependence on luminosity threshold to 
build catalogues for L > L* , L > 2L* , and L > 3L* galax- 
ies. Eor 2 = 1 catalogues, we used the HOD parameters mea- 
sured by Coupon et al. (201 1) and Zheng et al. (2007) respectively 
in the CFHTLS and DEEP2 surveys. In that case, while the two 
analysis use slightly different selection magnitude bands (g' and 
B), we assumed here that the later give comparable luminosities 
(the two bands largely overlap in wavelength) and interpolate be- 
tween the parameters indifferently of the band. The parameter Mq 
is poorly constrained by current observations and we decided to fix 
Mo = Mmin- We found that this approximation does not introduce 
any significant effect on the halo occupation and predicted cluster- 
ing. The characteristic absolute magnitudes at 2 = 0.1 and 2 = 1 
that we used to define the galaxy samples were taken from Blan- 
lon ct al. (2003) and llbcn ct al. (2005) respectively. The catalogue 
properties and HOD parameters are summarised in Table B 1 . 

This paper has been typeset from a TgX/ KlgX file prepared by the 
author 
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Table Bl. HOD parameters and associated galaxy catalogue properties. 



Redshift 


Luminosity 


Absolute magnitude 












Galaxy number 


2 


threshold 


threshold lo| 




•^log m 


log Mo 


log Ml 


a 


density 


0.1 


L> L* 


Mr - 51og(/i) < -20.44 


12.18 


0.21 


12.18 


13.31 


1.08 


0.347 


0.1 


L > 2L* 


Mr - 51og(/i) < -21.19 


12.89 


0.68 


12.89 


13.89 


1.17 


0.105 


0.1 


L > 3L* 


Mr - 5Iog(h) < -21.63 


13.48 


0.70 


13.48 


14.35 


1.23 


0.026 


1 


L> L* 


^B/g' - 51og(/i) < -20.78 


12.29 


0.35 


12.29 


13.36 


1.23 


0.200 


1 


L > 2L* 


Mg/g, - 51og(/i) < -21.53 


12.67 


0.38 


12.67 


13.77 


1.42 


0.072 


1 


L > 3L* 


Ms/g, -51og(h) < -21.97 


12.94 


0.40 


12.94 


14.10 


1.54 


0.032 


Masses are 


given in 


M0 and galaxy number densities in 10 


"2 Mpc 


-3 
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